Mon. Not. R. Astron. Soc. 000,[TH5](2008) Printed 27 July 2009 (MN BTgX style file v2.2) 



Quiet but still bright: XMM-Newton observations of the soft 
gamma-ray repeater SGR 0526-66* 

g : A. Tiengo/f P. Esposito, 12 S. Mereghetti, 1 G. L. Israel, 3 L. Stella, 3 R. Turolla, 45 
O S. Zane, 5 N. Rea, 6 D. Gotz, 7 and M. Feroci 8 

1 INAF/lstituto di Astrofisica Spaziale e Fisica Cosmica - Milano, via E. Bassini 15, 20133 Milano, Italy 
2 INFN - Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, via A. Bassi 6, 27100 Pavia, Italy 
3 INAF/Osservatorio Astronomico di Roma, via Frascati 33, 00040 Monteporzio Catone, Italy 
4 Universitd degli Studi di Padova, Dipartimento di Fisica, via F. Marzolo 8, 35131 Padova, Italy 
5 University College London, Mullard Space Science Laboratory, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK 
6 University of Amsterdam, Astronomical Institute Anton Pannekoek, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 
7 CEASaclay, DSM/lrfu/Service dAstrophysique, Orme des Merisiers, Bat. 709, 91191 Gif sur Yvette, France 
8 INAF/lstituto di Astrofisica Spaziale e Fisica Cosmica - Roma, via Fosso del Cavaliere 100, 00133 Roma, Italy 

Accepted 2009 July 25. Received 2009 July 24; in original form 2009 June 1 1 

ABSTRACT 

SGR 0526-66 was the first soft gamma-ray repeater (SGR) from which a giant flare was de- 
tected in March 1979, suggesting the existence of magnetars, i.e. neutron stars powered by the 
decay of their extremely strong magnetic field. Since then, very little information has been ob- 
tained on this object, mainly because it has been burst-inactive since 1983 and the study of its 
persistent X-ray emission has been hampered by its large distance and its location in a X-ray 
bright supernova remnant in the Large Magellanic Cloud. Here we report on a comprehensive 
analysis of all the available XMM-Newton observations of SGR 0526-66. In particular, thanks 
to a deep observation taken in 2007, we measured its pulsation period (P = 8.0544 ± 0.0002 
s) 6 years after its latest detection by Chandra. This allowed us to detect for the first time a 
significant reduction of its spin-down rate. From a comparison with two shorter XMM-Newton 
observations performed in 2000 and 2001, we found no significant changes in the spectrum, 
which is well modelled by an absorbed power-law with Nr = 4.6~^q 5 x 10 21 cm~ 2 and 
r=3.27t° o4. The high luminosity (~4x 10 35 erg s _1 , in the 1-10 keV energy band) still ob- 
served ^25 years after the latest detection of bursting activity places SGR 0526-66 in the 
group of bright and persistent magnetar candidates. 

Key words: ISM: individual: N49 - stars: neutron - supernova remnants - X-rays: individual: 
SGR 0526-66 - X-rays: stars. 



1 INTRODUCTION 

On 1979 March 5, an extremely bright gamma-ray burst (GRB), 
followed by a >60 s long tail puls ating at a period of 8.1 ± 0.1 s, 
was detected by many spacecrafts iMazets et alJI 19791) . The event 
was localized within the young (^5,000-10,000 years old) su- 
pernova remnant (SN R) LHA 120-N49 (N49) in the Large Mag- 
ellanic Cloud (LMC; Idine et alj 1 19821) . These properties indi- 
cated that the burst was emi tted by a young neutr on star, leading 
iDuncan & Thompson! J 19921) and Paczvnski J 1 9921) to propose the 
existence of neutron stars with magnetic fields of ~10 15 G, that 
were called magnetars. The detection of many weaker bursts from 

* Based on observations obtained with XMM-Newton, an ESA science mis- 
sion with instruments and contributions directly funded by ESA Member 
States and NASA. 

f E-mail: tiengo@iasf-milano.inaf.it 



the same direction in the following 4 years indicated that the March 
5 event was not a typical GRB, but an exceptional outburst from a 
small class of sources which had been just discovered and called 
soft gamma-ray repeaters (SGRs). Indeed the 1979 March 5 event 
from SGR 0526-66 was the first "giant flare" observed from a SGR. 
Only two other such events have bee n observed in the follow ing 
years, each one from a different SGR faurlev et alJll999Ll2005h . 

Up to now, only six SGRs have been discovered (plus a few 
candidates). They are characterised by the emission of short bursts 
of gamma-rays during sporadic periods of activity. In addition, they 
are also observed as pulsating X-ray sources with periods in the 2-9 
s range and persistent luminosities up to ~10 36 erg s" 1 . The mag- 
netar model was developed to explain both t heir bursting and per- 
sistent emission (Thompson & Duncan! 1 19951) and later extended 
(Thompson & Duncan 1996) to the interpretation of the anomalous 
X-ray pulsars (AXPs). These are a group of 9 X-ray sources (plus 
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some candidates) with similar properties to those of the SGRs (see 
lMereghettill2008l for a review). 

The persistent X- ray emission from SG R 0526-66 was first 
detected with ROSAT dRofhschild et al.lll994l) and then observed 
by C handra in 2000 and 2001 jKulkarni et all 120031 : IPark et all 



Table 1. Log of the XMM-Newton observations of SGR 0526-66. 



2003), with a constant X-ray luminosity of ~10 36 erg s 1 (un- 
absorbed, in the 0.5-10 keV energy range). Being SGR 0526-66 
a rather faint X-ray source, the pulsation of its persistent emis- 
sion was detected on ly in the two Chand ra observations carried 
out in 2000 and 2001 (Kulkarni et al. 2003). The periods measured 
with Chandra are compatible with the 8 s period detected in the 
pulsating tail of the 1979 giant flare and give a spin-down rate of 
P — (6.5 ± 0.5) x 10 _1 s s~ , a value in the same range as that 
of the other magnetar candidates. 

In order to measure again the pulsation period and search for 
long-term flux and spectral changes, we obtained a ~70 ks long 
observation of SGR 0526-66 with XMM-Newton, 6 years after the 
latest X-ray observation. In the following we report on the results 
of this recent observation, together with the analysis of two short 
archival XMM-Newton observations performed in 2000 and 2001. 



Obs. 


Date 


Instrument 


Mode a 


Filter 


Net exposure 


A b 


2000-07-08 


PN 


FF 


medium 


4.3 ks 






MOS1 


FF 


medium 


8.4 ks 






MOS2 


FF 


medium 


8.5 ks 


B 


2001-04-08 


PN 


FF 


medium 


6.7 ks 






MOS1 


FF 


thick 


4.6 ks 






MOS2 


FF 


medium 


12.6 ks 






PN 


SW 


medium 


5.8 ks 






MOS1 


FF 


medium 


12.1 ks 


C 


2007-11-10 


PN 


LW 


thick 


60.3 ks 






MOS1 


SW 


thick 


69.7 ks 






MOS2 


SW 


thick 


69.7 ks 



a The time resolution of the operating modes are: PN Full Frame (FF): 73 
ms; PN Large Window (LW): 48 ms; PN Small Window (SW): 6 ms; MOS 
Full Frame (FF): 2.6 s; MOS Small Window (SW): 0.3 s. 
b 6' off-axis. 




2 OBSERVATIONS AND DATA ANALYSIS 

SGR 0526-66 was observed by XMM-Newton on 2007 Novem- 
ber 10 for about 70 ks. The field containing SGR 0526-66 had 
already been observed by XMM-Newton, with shorter exposure 
times, on 2000 July 8 and on 2001 April 8. In the 2000 observa- 
tion SGR 0526-66 was ~6' off-axis, while in the other observa- 
tions it was on-axis. We concentrate here on the analysis of the 
data collected with the EPIC instrument, which is co mposed by a 
PN dStriider et alj200lh and two MOS X-ray cameras dTurner et aT] 
1200 C sensitive in the 0.2-15 keV energy range. Details on the in- 
strument settings (optical blocking filter and operating mode) for 
each observation are listed in Table Q] For the longest observation 
we used al so the data collected by the Reflection Grating Spectrom- 
eter (RGS. Iden Herder et alj200lh . which worked in parallel to the 
EPIC instrument and had a net exposure time of 7 1 ks for each of 
its two units (RGS1 and RGS2). This high resolution spectrometer 
is sensitive in the 0.35-2.5 keV energy range. 

All the data were processed using the XMM-Newton Science 
Analysis Software (SAS version 8.0.0) and the calibration files re- 
leased in August 2007. The standard pattern selection criteria for 
the EPIC X-ray events (patterns 0-4 for PN and 0-12 for MOS) 
were adopted. The RGS analysis followed the standard selection 
criteria as wellQ Response matrices and ancillary files for each 
spectrum were produced using the SAS software package and the 
spectra were fitted using XSPEC version 1 1.3.1. All errors reported 
in the following analysis are at 1 a. 




Figure 1. EPIC PN image of the field of SGR 0526-66 during the 2007 
observation. Photon energy is colour-coded: red corresponds to 0.5-2 keV, 
green to 2—4 keV, and blue to 4-10 keV. The source and background regions 
considered for the analysis are overplotted. 



in the fits with a free normalization and a fixed spectral shape de- 
termined as explained below. To reduce the contamination from the 
soft X-ray emission of N49, we extracted the SGR EPIC spectrum 
in the 1-10 keV energy range from a circle of 10" radius (this in- 
cludes 60% of the point source counts at 5 keV). The background 
spectrum was extracted from a region outside the SNR, but in the 
same CCD as the SGR (see Figure[T|l. 

To model the soft and line-rich spectrum of the SNR, we 
took advantage of the high-resolution spectra collected by the 
RGS instrument during the longest observation. We extracted the 
first order spectra from the standard region normally used for 
point sources, setting the centre of the SNR as source position 
(such a selection includes most of the photons detected from the 
http://xmm.esac.esa.int/external/xmm_user_support/documentation/sas_usg/USG/ SNR, thanks to its relatively small extension). The RGS spectral 



2.1 Spectral analysis 

The spectral analysis of SGR 0526-66 with XMM-Newton is com- 
plicated by the location of this source within the bright SNR N49, 
whose spatial extent (~40" radius) is only slightly larger than the 
instrumental point-spread function (the 90% encircled energy frac- 
tion for a point source is ~40"). Rather than attempting to subtract 
the SNR contribution as a background component, we included it 



See 
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analysis was restricted to the 1-2 keV energy range. To model 
the SNR above 2 keV, we extracted a 1-10 keV PN spectrum 
from a 40" circle centred in the middle of the SNR and fitted 
it together with the spectra of the two RGS un its. Based on the 
results of previous X-ray observations of N49 JParketalj|2003l ; 
lBilikovaetaT]|2007h . we used a model consisting of the sum of 
two plane-parallel shock components at different temperatures 
(VPSHOCK in XSPEC) both corrected for photoelectric absorption 
(PHABS in XSPEC). To this we added an absorbed power-law 
to account for the emission from SGR 0526-66. Its parameters 
were fixed at the best-fit values found with Chandra : Nu 
=5.6xl0 2 1 cm" 2 , r=3.06, no rm=l. 18x 10~ 3 s" 1 cm" 2 keV" 1 
at 1 keV I Kulk arni et al.l I2003T) . An overall normalization factor 
for each spectrum was also included to account for the cross- 
calibration uncertainties between the two RGS and the EPIC PN. A 
satisfactory fit (x 2 =2170/935 de grees of freedom, d.o.f. ]J was ob- 
tained with the following parameters: Ah =(1.3±0.3)x 10 21 
cm -2 , fcTi=0.577t°°° 2 keV, ri=5.4±J;£xl0 12 s cm" 3 , 
fcT 2 =1.10±0.01 keV, t 2 >3x10 13 s cm" 3 , Ne/Ne Q =0.66±0.02, 
Mg/Mg Q =0.59±0.01, Si/Si e =0.79±0.01, S/S Q =1.00±0.04, 
Ca/Ca Q =0.6±0.4, Fe/Fe Q =0.45±0.02. The abundances of the 
other elements are fixed to the Solar values dAnders & Grevessel 
1989). These parameter s are in good ag r eement with previou s 
studies of this SNR jPark et all l20oj iBilikova et ail l2007h . 
The unabsorbed flux in the 1-10 keV energy is 7.1xl0~ 12 
erg cm" 2 s" 1 for the SNR and 1.6xl0" 12 erg cm" 2 s" 1 for the 
SGR. 

We then fitted the EPIC spectra extracted from the small re- 
gion around SGR 0526-66 with an absorbed power- law3 plus the 
model of the SNR described above. All the SNR model parame- 
ters were fixed at their best-fit values, except for the normaliza- 
tion, in order to properly account for the unknown intensity of 
the SNR emission in the source extraction region0 A good fit 
(X 2 = 274.1/233 d.o.f., see Figure^ is obtained with a hydro- 
gen column density Ah = (4.6^0 5) xlO 21 cm" 2 and a photon 
index r=3.27+{J'°,4. The lack of systematic residuals in correspon- 
dence with the SNR brightest spectral lines (see Figure O indi- 
cates that the SNR contamination is sufficiently well modelled. 
The 1-10 keV (unabsorbed) flux of the power-law component is 
(1.25^q qj) x 10~ 12 erg cm" 2 s" 1 . This corresponds to a luminos- 
ity of 4.3 x 10 35 erg s" 1 for a distance of 55 kpc. 

Since no time variability is expected in the SNR contribution 
we can fit the SGR spectra of the older EPIC observations with 
the model described above, keeping the SNR model normalization 
fixed at the value obtained in the longest observation. The best-fit 
parameters for an absorbed power-law model are reported in Ta- 
ble^ No spectral variability is detected and significant (>3<r) flux 
variations larger than ~50% among the different XMM-Newton ob- 
servations can be excluded. 



2 Although this fit is not statistically acceptable, a detailed modelling of 
the SNR emission is beyond the scope of this paper and so we did not adopt 
more complex spectral models. 

3 More complex spectral models, which are usually used to fit magnetar 
spectra, were not adopted in this case due to the uncertainty of the back- 
ground subtraction. 

4 We also applied a second normalization factor to both the SNR and the 
SGR models to account for the cross-calibration uncertainties between the 
EPIC cameras; the maximum flux discrepancy we find between the EPIC 
cameras is lower than 15%. 
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Figure 2. EPIC PN spectrum of SGR 0526-66 collected during the longest 
XMM-Newton observation in 2007. The model, obtained by a simultaneous 
fit of the PN and MOS spectra (see the corresponding parameters in Ta- 
ble [2), is an absorbed power-law (green) and the sum of two plane-parallel 
shock components (in blue the warmer and in red the cooler one) to model 
the contamination from the SNR. 

2.2 Timing analysis 

For the timing analysis the photon arrival times were converted to 
the solar system barycenter with the BARYCEN SAS task. We used 
the same extraction region adopted for the spectra, but we per- 
formed the analysis in the 0.65-12 keV energy band to optimize 
the signal to noise ratio. The pulse periods measured with Chandra 
using the Z% test were 8.0436(2) s on 2000 l anuary 4 and 8.0470(2) 
s on 2001 August 31 iKulkarni et al1l2003ri . We searched for pul- 
sations in the period range 8.0464-8.2423 s, selected by consid- 
ering the 3a lower limit on the most recent spin period value re- 
ported by Kulka rni et alj d2003h and extrapolating to the epoch of 
the XMM-Newton observation under the assumption of a (conser- 
vative) period derivative of ^ P ^ I0" 9 s s" 1 . In Figure[3]we 
show the Z\ periodogram obtained by using the combined events 
from the PN and MOS cameras. Taking into account the number of 
searched periods (419), the peak value of 31.08 (chance probabil- 
ity of ~3 x I0" 6 ) at ~8.0544 s corresponds to a significance of 
~99.88% (that is a 3.2<r detection). 

To better estimate the period, we used an epoch folding tech- 
nique and fitte d the peak in t he \ 2 versus trial period distribution 
as described in lLeahvl dl987l) . This led to P = 8.0544 ± 0.0002 
s. The resulting EPIC folded lightcurve is shown in Figure [3] The 
pulsed fraction, defined as (C ma x — C m i n ) / (C max + C m i n ), where 
C max and C m in are the background-subtracted count rates at the 
peak and at the minimum, is (13.6 ± 0.9)% in the 0.65-12 keV 
energy range 0We repeated a similar analysis on the data-sets from 
the two short XMM-Newton observations carried out in 2000 and 
2001 (see Table [T), but, due to the paucity of photons, we did not 
detect any significant periodicity. 

Considering also the two periods measured by Chandra 
jKulkarni et al.l [2003), a linear fit to the period evolution of 
SGR 0526-66 is unacceptable { X 2 of 15.6 for 1 de gree of free- 
dom). On the other hand, the period derivative inferred from our 
period measurement and the most recent Chandra detection is 
(3.8 ± 0.1) x 10" 11 s s~ . Comparison with the value derived 

5 The SNR contamination, evaluated from the model described in Section 
2.1, is also included in the background. The error in the pulsed fraction does 
not include the systematic uncertainty due to the SNR contamination. 
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Table 2. Best-fit spectral parameters for the three EPIC observations (see Table[T} of SGR 0526-66 in the 1— 
10 keV energy range. The spectral model consists of a fixed component modelling the SNR contamination 
(see text for details) and an absorbed power-law model for the SGR emission. 



Observation 


SNR norm. a 


Nh 
(10 21 cm" 2 ) 


r 


Flux b 

(10 -12 erg cm -2 s _1 ) 


Xl (d-o.f.) 


A 
B 
C 


0.117 (fixed) 
0.117 (fixed) 
1 17 +0.002 
uai '-0.005 


J - J -l.l 

4 6+0.7 


3 l +0 - 4 

1 1+0.1 

J.o_ 2 
i 97 +0.07 
°- z '-0.04 


i i+0.2 

L -°-o.i 
1.3+0.1 

, 9 <r+0.05 
l -"-0.02 


1.04(65) 
0.88 (166) 
1.18(233) 



a Normalization factor applied to the best-fit model of the RGS and PN spectrum of the whole SNR. 
b Unabsorbed flux in the 1-10 keV range. 




sphere and should be associated to a low rate of bursting activity, a 
spe ctral softening and a de crease of the persistent X-ray luminos- 
ity dThompson et alj|2002l) . This behavior is sometime s observed 
in magnetar candidates (see, e.g.. lMereghetti et"al]|2005l), but s ome 



Gavriil & Kaspll2004l) . In 



Figure 3. Z?, diagram of the long XMM-Newton observation (PN and MOS 
data in the 0.65-12 keV energy range) of SGR 0526-66 in the range used 
for the period search (see Section 2.2). The peak at 8.0544 s is significant at 
3.2cr. Inset: The corresponding EPIC pulse profile (0.65-12 keV, not back- 
ground subtracted). 



from the two Chandra measurements (P — (6.5 ± 0.5) x 10 11 s 
s _1 ) indicates a significant decrease of the spin-down rate. 



3 DISCUSSION AND CONCLUSIONS 

Due to its location in the LMC, SGR 0526-66 is one of the mag- 
netars less frequently observed. In addition to its larger distance 
with respect to Galactic magnetars, the analysis of its persistent X- 
ray emission is complicated by the bright soft X-ray emission of 
the N49 SNR. With a deep XMM-Newton observation performed 
in 2007, we could measure the pulsation period of SGR 0526-66 
which was previously detect ed only in the pulsa ting tail following 
the 1979 March 5 giant flare dMazets et aljl979h and in two Chan - 
dra observations taken in 2000 and 2001 dKulkarni et alj|2003h . 
The pulsation profile (shown in Figure^ is double-peaked and the 
pulsed fraction is (13.6±0.9)%. Although, to our knowledge, the 
pulse profile of SGR 0526-66 was never published before, these 
might be permanent properties of this source since a non-sinusoidal 
modulation and a pul sed fraction around 10% were also reported 
for the Chandra data dKulkarni et alj2003l) . 

Although the period measurements of SGR 0526-66 are very 
sparse, the value we derived shows, for the first time, a signifi- 
cant decrease in the spin-down rate of this source. In the mag- 
netar model, a reduction of the spin-down rate can be interpreted 
as an indication of a more relaxed state of the twisted magneto- 



exceptions have been found (see, e.g. 
the case of SGR 0526-66, the bursting activity is indeed very low, 
since no bursts have been detected since 1983. However, we note 
that some bursts might have been missed due to its large distance 
and its location in a sky region not frequently monitored by 7-ray 
instruments. 

The X-ray luminosity measured in 2007 is ~ 30% lower (and 
the spectrum slightly softer) than that repo rted from the analysi s 
of Chandra data taken in 2000 and 2001 dKulkarni et ai]|2003b . 
However, due to the different characteristics of the instruments and 
additional uncertainties due to the presence of contamination from 
the SNR diffuse emission, we consider these changes well within 
the systematic uncertainties. Using the two shorter XMM-Newton 
observations, taken almost simultaneously with the Chandra ones, 
we can instead extract the spectrum and model the SNR emission 
in the same way as we did for the 2007 observation. In this case we 
are dominated by statistical errors and no significant changes in the 
spectral shape and source flux are detected. 

The luminosity of SGR 0526-66, which can be well deter- 
mined thanks to its accurately know n distance, is higher than that o f 
most magnetar candidates (see e.g. lDurant & van Kerkwijk 2006). 
Since this high luminosity has not substantially varied for at least 
several years, it is probably a long-lasting property of this magnetar 
rather than a transient bright state related to its past bursting activ- 
ity, culminating with the giant flare of 1979 March 5. This behavior 
is radically different from the one displayed, as an example, by 
SGR 1627-41 that, after two distinct periods of bursting activity, 
rapidly decreased its persiste nt X-ray luminosity down to ~10 33 
erg s _1 dEsposito et alj|2008t) . While it is unclear whether this be- 
havior is related to intrinsic differences between the sources (e.g. 
the magnetic field) or different evolutionary stages, this seems to 
support the emerging trend of separating the magnetars into tran- 
sient and persistent objects, rather than in AXPs and SGRs. In this 
framework, SGR 0526-66 should therefore be considered a mem- 
ber of the persistent magnetars group. 
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